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1 Introduction 



In the studies that are now being done to prepare for physics at the LHC, many 
new approaches have been proposed to the old problem of generating parton showers. 
The workhorse event generators PYTHIA [I] and HERWIG [2] generate parton show- 
ers by successive radiations from individual partons. The 'splitting functions' that 
define the radiation pattern are taken to be the kernels in the Altarelli-Parisi equa- 
tion [3_,4J. This guarantees that the radiation pattern is correct in the region in which 
two partons become coUinear. Marchesini and Webber pointed out that it is also 
important to include color interference between emissions from different partons [5]. 
In the workhorse generators, this is implemented by angular ordering of emissions. 

The program ARIADNE, by Andersson, Gustafson, Lonnblad, and Pettersson, 
took a different approach, implementing color coherence by considering the QCD 
dipole to be the basic object that radiates a parton [Gj^. The basic branching pro- 
cess in a parton shower is then a splitting in which two partons forming a color dipole 
radiate a third parton. This approach has been taken up recently by a number of 
authors. It is the basis for the VINCIA shower by Giele, Kosower, and Skands [8] and 
the parton shower implementation in SHERPA by Krauss and Winter [H]. We are 
also developing a parton shower based on this approach [TU]. In the years between 
ARIADNE and the newer works, the term 'dipole' has been applied in QCD to a 
different strategy based on 1 — )■ 2 splittings with recoil taken up by a third parti- 
cle [11]. To avoid confusion, we will follow |8] in calling the initial two-parton state 
an 'antenna' and a branching process with 2 — )■ 3 splittings an 'antenna shower'. 

Central to the antenna shower is the 2 — )■ 3 splitting function, the function that 
gives the relative branching probabilities as a function of the final momenta. The 
original ARIADNE program used an ad hoc proposal satisfying the basic consistency 
requirements. It would be better to have a prescription that can be directly derived 
from QCD. Splitting to three partons has been studied in great detail in the QCD 
literature, but not for this application. Collinear systems of three partons are a part 
of the infrared structure of QCD at next-to-next-to-leading order, and calculations 
that reach this level need an explicit prescription for treating this set of infrared sin- 
gularities. Kosower [12] defined the 'antenna function' as a basic starting point for the 
analysis of this problem. Many authors have computed antenna functions [T3lll4|15] . 
Quite recently, Gehrmann-De Ridder, Gehrmann, and Glover have built a complete 
formalism of 'antenna subtraction' for NNLO calculations [16]. The kernel in their 
theory can be interpreted as a 2 — )• 3 splitting function, and it has been used to 
perform 2 — )■ 3 splitting in the VINCIA shower [8]. 

In this paper, we will take a much more direct route to the construction of 2 — )■ 3 
splitting functions. We will compute these functions by writing local operators that 
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create two-parton final states and computing their 3-parton matrix elements. These 
calculations are very straightforward. They can be used to treat individually all 
possible sets of polarized initial and final partons. 

This paper is organized as follows: In Section 2, we will present our complete set 
of spin-dependent 2 — )■ 3 splitting functions. In Section 3, we will give the derivation 
for the cases with total spin zero. In Sections 4 and 5, we will give the derivation for 
the cases with nonzero total spin. 

All of these derivations will be done in the kinematics of final-state radiation. 
This is the easiest situation to visualize and understand. However, the same splitting 
functions can be used, after crossing, to describe parton emissions that involve initial- 
state particles. We will explain how to use our expressions for initial-state showers in 
Section 6. 

The 1 — )■ 2 Altarelli-Parisi splitting functions are universal in the sense that they 
result from a well-defined singular limit of QCD amplitudes. For 2 — )■ 3 splitting func- 
tions there is no such universality. The coUinear and soft limits must agree with the 
known universal values, but away from these limits there is no unique prescription. 
Earlier in this introduction, we made reference to a number of previous proposals for 
the spin-averaged antenna splitting functions. All of these, including the ARIADNE 
splitting functions, have the correct soft and collinear limits and so satisfy the basic 
requirements. In Section 7, we will give a detailed comparison of the 2 — )■ 3 split- 
ting functions obtained using our method to previous proposals for these splitting 
functions. 



2 Proposal for the 2^3 splitting functions 

We begin by defining variables for 2 — 3 splitting. There are three cases of 
splittings that are needed for antenna showers: the final-final (FF) splitting, in which 
a third particle is created by coherent radiation from a two-particle system in the final 
state; the initial-final (IF) splitting, in which a third particle is created by coherent 
radiation from an initial- and a final-state particle; and initial-initial (II) splitting, in 
which a third particle is created by coherent radiation from two initial-state particles. 
It is easiest to understand the kinematics of antenna splitting for the FF case. In 
this section, we will explain this kinematics and give a precise prescription for the 
splitting functions. In Section 6, we will extend our prescription to the IF and II 
cases, in such a way that the same splitting functions can be used in those cases. 

Consider, then, a two-parton final-state system {A, B) that splits to a 3-parton 
system (a, c, 6), conserving momentum, as shown in Fig. fl](a). Let Sij = [ki + /c^)^. 
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Figure 1: (a) Kinematics of 2 — )• 3 splitting in the final state (FF) case, (b) Phase space 
for 2 — )• 3 splitting in the FF case. The six regions corresponding to different orderings of 
■Sab) Saci She are shown. The region that should be well described by an antenna splitting 
AB — > ach is shaded. 



and let Q = kA + ks = ka + + k^. 

The fractional invariant masses in the final state are 

Sab Sac Sfjc , , 

Vab = , Vac = , Vbc = ■ (1) 

Sab Sab sab 

The momentum fractions of the three particles in the (AB) frame are 

2Q ■ ka 2Q ■ kh 2Q ■ k^ , < 

) ^b ) ■ ) 

Sab Sab sab 

These obey the identities 

Vab = (1 - Zc) , Vac = (1 - Zb) , Vbc = (1 " Za) ■ (3) 

and 

Vab + Z/ac + Z/6c = 1 , Za + + = 2 . (4) 

The FF phase space covers the triangle Za < 1, Zb < 1, Za + Zb > 1. We can 
divide this phase space into six triangles, each of which has a different ordering of 
the three quantities yab, Vac, Vbc, as shown in Fig. [T|b). An antenna shower should 
give an accurate description of the dynamics in the two regions i/ac < Vbc < Vab, 
Vbc < Vac < Vab that are shaded in the figure. 

A general problem in the generation of QCD radiation is that of possible double- 
counting. Consider, for example, the process e~^e~ — t- qggq. In some part of the phase 
space, the first g can be considered to be radiated from the antenna of the q and the 
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second g; in another, the second g can be considered to be radiated from the first 
g and the q. These regions should be disjoint in the full 4-body phase space. The 
complete solution to the problem is beyond the scope of this paper. In simple terms, 
though, we can make the separation by choosing the radiated gluon to be softer than 
the gluon from which it radiates. This corresponds to integrating each antenna only 
over the shaded region in Fig. [l|b). A similar approximate solution to the double- 
counting problem will apply in the other kinematic regions discussed in Section 6. A 
more detailed discussion of this issue can be found in PITU] . 

Radiation from different QCD antenna is strictly independent and non-interfering 
only in the limit of a large number of colors in QCD, Nc ^ 1. Keeping only terms 
leading in Nc is known to be a good approximation to full QCD in many circum- 
stances. In particular, parton shower algorithms are correct only to leading order in 
Nc- In this paper, we will explicitly work only to the leading order for large N^. 

In the limit of large the rate for a 2 — > 3 splitting is given by a formula of the 
form 

Nc'^ J dZadZb ■S{Za,Zb,Zc) (5) 

For example, in e~^e^ -t- q-g+q^, 

do dzadzh "An (1 - - Zb) 

where (a, c, 6) are the {q,g,q), respectively, — and + denote left- and right-handed 
helicity, and ctq is the cross section for e+e~ — t- [I7j. Eq. ([5| will be our basic 

formula of reference. Using this notation, we can write the various 2 — )■ 3 splitting 
functions as 

^ ^ AfjZq, Zb, Zc) 

y ahV acVhc 

where the numerator is a simple function of the Zi. For example, for the splitting 
— )• given above, 

N = Vabzl = (1 - Z,)zl . (8) 



In Table [T| we give our proposal for the numerator functions for all possible cases 
of massless quark and gluon splittings. The expressions are all monomials in the Uij 
and Zj. 

In the FF kinematics, all of the yij and Zi are positive and so S{za, Zb, Zc) in ([T]), 
is always positive. In IF and II kinematics, some Uij and Zi will become negative. In 
most cases, the correct prescription is to take S{za, Zb, Zc) to be the absolute value of 
the expression in Table [1} However, there is a line within the IF region where Za or 
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Table 1: Numerator functions M{za^ Zf,, Zc) for the spin-dependent 2 — )• 3 splitting functions 
AB — )• ach: S = M / {yabVacVbc)- Each line gives a choice of AB. The labels denote the 
polarization of the three final particles with the radiated particle c in the center: (/iq, /ic, /ife). 
The empty columns are forbidden by quark chiral symmetry. By the P and C invariance of 
QCD, the same expressions apply after exchanging — <-)•+, g o or ABacb -f-)- BAbca. 



Zh crosses from positive to negative values. A few entries in the Table change sign 
across this line. We recommend that those entries be set to zero when Za or Zfy are 
negative. We will give a detailed discussion of these points in Section 6. 

The splitting functions S must give the correct universal behavior in the soft and 
collinear limits. In the soft limit, Zc — ?• 0, the numerators must go to 1 if the flavor 
and helicity of the final partons a and b match those of the initial partons A and B; 
otherwise, the numerators must go to 0. It is easy to check that this test is satisfied. 

In the collinear limits, we will insist that each antenna has the collinear behavior 
required in QCD. One often hears the following statement about soft and collinear 
limits: In dipole splitting (1 — )■ 2 emission), each dipole has the correct collinear 
behavior but the correct soft behavior is obtained by combining neighboring dipoles. 
In antenna splitting (2 — )■ 3 emission), each antenna has the correct soft limit but the 
correct collinear behavior is obtained by combining neighboring antennae. However, 
in the large Nc limit, which we take to guide our intuition, different antennae are 
independent radiators with different, non- interfering, colors flowing in them. From the 
viewpoint of this limit, each antenna, separately, must give both the correct pattern 
of soft radiation and the correct pattern of collinear radiation. This philosophy differs 
from that of the ARIADNE group pi7] and of [H]. We will discuss this point further 
when we compare with their results in Section 7. 

The collinear radiation from a given hard gluon is then the sum of two contri- 
butions, one from each of the two antennae to which that hard gluon belongs. In 
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Table 2: Spin-dependent Altarelli-Parisi splitting functions P{z) for splittings B — )• cb. 
The labels denote the polarization of the two final particles with the radiated particle first: 
{hc,hh). The empty columns are forbidden by quark chiral symmetry. By the P and C 
invariance of QCD, the same expressions apply after exchanging — o + or (7 -f-)- g. 



the large limit, these correspond to radiation from the color and anticolor lines of 
the gluon. A single antenna, which has one of these contributions, then has ^ of the 
standard coUinear emission rate. This factor of | enters the check will we perform 
in a moment. The factor comes entirely from bookkeeping and is independent of the 
question of double-counting discussed briefly earlier in this Section. 

We now discuss the check of coUinear limits. Consider the limit in which c becomes 
coUinear with a. In this limit, 

Zc^ Z , Za ^ {1 - Z) , Zb^ 1 , Vac^O . (9) 

The 2 — 7- 3 splitting function must reduce to 

S ^ —P{z) , (10) 

Vac 

here P{z) is the relevant spin-dependent Altarelli-Parisi splitting function. These 
were presented in the original Altarelli-Parisi paper [3] and are reviewed in Table [2] 
The functions are normalized as in ([s]), and as described in the previous paragraph: 
We take the large Nc limit and divide by 2 where necessary to give the contribution 
from one QCD antenna. The denominator of ([T]) tends to yacz{l — z) in this limit. 
Then it is easy to check that the numerators match correctly in all cases. The limit 
in which c becomes coUinear with h can be checked in the same way. 

When the coUinear limits and the soft limit are all nonzero, there is a unique 
monomial of the y's and z's that gives all limits correctly. In the other cases, there is 
some ambiguity. In all cases, it would be desirable if the results in Table [T] could be 
derived directly by simple Feynman diagram computations. In the next few sections, 
we will present those derivations. 
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3 Spin-0 case 



To compute the 2 — )■ 3 splitting functions, we will use the following method: Write 
an operator that, at the leading order, creates a 2-parton state with definite helicity. 
Then, compute the 3-particle matrix element. This realizes in a very simple way the 
splitting process illustrated in Fig. [Tj 

To create massless quarks and antiquarks of definite helicity, we will use the ap- 
propriate chiral fermion fields. To create gluons of definite helicity, we will use the 
operators 

a-F= ^cr"^cf"F^„ , a-F= ^a'"a"F^„ , (11) 

where cr"*, are the 2x2 matrix entries of the Dirac matrices in a chiral basis and 
Fmn is the gluon field strength tensor. At leading order, a ■ F creates a + helicity 
gluon, and a ■ F creates a — helicity gluon. 

The 2-parton state g+g+ in the first line of Table [l] can be created from the spin-0 
operator 

O = hrlia . Ff] . (12) 

We can then compute the splitting function for this polarized initial state explicitly 
from the definition 



M{0^ ach) ^ 



(13) 



M{0^ AB) 

In the next few sections, we will compute all of the splitting functions in Table [T] 
using this formula, with a different choice of the operator O for each line of the table. 



To evaluate (13), we need to compute the matrix elements of (9, with total mo- 
mentum Q injected, to 3-gluon final states. The result can be expressed in terms 
of color-ordered amplitudes. We identify the color-ordered amplitude that multiplies 
the color structure tr [T'^T'^T^] with the splitting function. To carry out these compu- 
tations, we will use the spinor product formalism. That is, instead of working with 
4-vectors, we will use as our basic objects the spinor products 

(u) =u-{i)u+{j) , [ij] =u+{i)u^{j) . (14) 

These objects obey 

m? = m? = s^,. (15) 

Methods for QCD computations with spinor products and color-ordering are ex- 
plained in [T^PU] . In this notation, the matrix element for O to create a g+g+ final 
state is 

{g+g+\0\Q) = [ABf . (16) 
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Figure 2: Feynman diagrams for the computation of the gg — )• ggg spHtting functions. 



The three-gluon matrix elements of the operator ( 12 ) are given by the diagrams in 



Fig. [2j These diagrams have already been analyzed by Dixon, Glover, and Khoze as 
a part of their analysis of the coupling of the Higgs boson to multi-gluon states [20] . 
They find 



AiO ^ g+g+g+] 
A{0^g+g+g-] 
AiO ^ g^g-g^] 
A{0 -> g-g+g+: 



{ac){cb){ba) 

[ac] [cb] [ba] 

[ab]^ 
[ac] [cb] [ba] 

[ac] [cb] [ba] 



(17) 



and zero for the other four cases. After squaring, using (15), and dividing by the 
square of (16), we obtain the first line of Table [l] 



One of the major points of |20] is that the results ( |17[ ) belong to series of Maximally 
Helicity Violating (MHV) amplitudes that have a simple form for any number of 
gluons emitted. Actually, all of the amplitudes that we will compute in this paper 
are similarly simple and belong to MHV series. The use of MHV amplitudes to study 
antenna splitting is explored for higher-order processes in [T5] . 

In principle, the initial state g+g+ could also have been created by an operator of 
spin 2, or some higher spin. This would have led to a more complicated expression 
for the 2 — )■ 3 splitting function, with, however, the same soft and coUinear limits. 
This illustrates the ambiguity in the definitions of 2 — )■ 3 splitting functions refered to 
in the introduction. The simplest results are obtained using the operator of minimal 
spin, and we will make that choice in all of the examples to follow. 



The diagram shown in Fig. [3] gives the splitting of the two-gluon initial state to 
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Figure 3: Feynman diagram for the computation of the gg — )• qqg spHtting functions. 



qqg. We find 

[abf 



A{0 q+q^g^ 



ac\ 

2 



AiO ^ q_q^g^) = (18) 
[ac\ 

There is no splitting to a final g^. This gives the result in the third line of the table. 
The initial state can also be created by a spin operator 

O = (iLqR . (19) 
The matrix element for this operator to create a q~q_ final state is 

O\0) = {AB). (20) 



{abf 



A straightforward calculation gives 



AiO ^ q-g+q_ 



(ac) {cb) 

Sab 
[ac] [cb] 

These give the results shown in the sixth line of the table 



AiO ^ q.g.q_) = (21) 



4 Spin-1 and spin-2 case 

In JE\ , the 2 — )■ 3 splitting function for qq — qgq was derived from the cross section 
for e'^e~ — )■ qgq. From the point of view of the previous section, this corresponds to 
creating the 2- and 3-parton final states using the operator 

O = q^^^qL . (22) 
9 



To obtain a definite matrix element, we must contract this operator with a polariza- 
tion vector. A convenient choice is to introduce two new massless vectors 1 and 2, 
such that ki + k2 = k^ + kB, and to choose the polarization vector to be = (1|7'^|2]. 
This is effectively the procedure of decaying the massive vector that couples to the op- 



erator (22) into a pair of massless vectors to facilitate the analysis; this is a standard 
method in spinor product calculations [21]. We then recast 



(23) 



O = t;Ql1"'Ql {l\lm\2] 



The matrix element of (|23) to a q^q_^_ state is 

O\0) = -{1A)[2B] . 



(24) 



The direction of the 1-2 system chooses the helicity of the final partons. In this case, 
there is only one choice, and so the amplitude vanishes when 1 is parallel to A or 2 is 
parallel to B. This will not always be true in our later examples. But, we will always 
be able to choose the desired helicity of A and B by choosing 1 parallel to B and 2 
parallel to A. 



The matrix elements for the operator (23|) to create 3-parton final states are 

AiO q-g+q+\ 

A{0 ^ g_(7_g+) = 7\:;7 . (25) 



(ac) (cb) 

[26]^(12) 
[ac] [cb] 



To compute the results in the fifth line of the table, we must essentially divide 



(25) by (24) and square the result. To do this, we need a prescription for treating the 
expressions (la) and [26] in the numerators. The problem of relating the vectors a, 
b, c to A and B in an antenna splitting was discussed at length by Kosower in [22] : 
that paper gives a general treatment in terms of reconstruction functions to provide 
expressions that can be smoothly integrated in higher-order QCD calculations. This 
discussion is generalized to the initial-state channels in [23]. Here, we will take a more 
ad hoc approach that leads to the simplest formulae with correct singular limits. 

Formulae for (la) and [2b] that are simple and become exact in the coUinear and 
soft limits are found by approximating a coUinear with A and b collinear with B. 
Then identifying 1 with B and 2 with A gives 

I (la) |2 = SBa ^ ZaSAB , \ (16) |' ^ , | (2a) 1^ ^0,1 (26) I' = SAb ^ ZbSAB , (26) 
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and similarly for the conjugate products. Using this prescription, one obtains the 
fifth line of the table. This is a more formal version of the argument for these entries 
already given in Section 2. 

In our calculations, we will encounter two more numerator objects that require 
reconstruction, namely, (Ic) and (2c). The prescription above gives 

\{lc)\^ = Sbc^ {ybc/Zh)sAB , \{'^C)\'^ = SAc^ {yac/Za)sAB ■ (27) 

However, it is potentially dangerous to write factors of Za-, zi, in the denominator. 
We will see in Section 6 that such factors would create unphysical singularities when 
continued to the IF kinematics. These unphysical singularities are avoided in the 
general formalism used in |22] , but at the price of introducing much more complicated 
formulae. Fortunately, we will see that (Ic) arises only in situations where there is 
no collinear singularity with c parallel to h. In such cases, the remaining universal 
singular terms — the collinear singularity with c parallel to a and the soft singularity — 
correspond to kinematic limits with zj, — 1. A similar consideration applies to (2c). 



Thus, we choose, instead of using (27), to evaluate these quantities as 



\{lc)\^ = S Be ^ VbcSAB , \{2c)\^ = SAc ^ VacSAB ■ (28) 

This gives an incorrect shape in a region where a and h are collinear, but, hopefully, 
we will not use the AB acb splitting function to evaluate the rate to fill this region 
of phase space. 

Another choice for evaluating (Ic) and (2c) is to replace both expressions by Zc- 
However, the spinor product (Ic) vanishes in the be collinear limit but not in the 
ac collinear limit, and conversely for (2c), so this choice does not give the universal 
singularities correctly. 

We now apply this formalism to compute the second and fourth lines of Table [T| 
associated with the g-g+ antenna. This antenna is created by the spin-2 operator 
tr[7"^(cf ■ -F)7"((T • F)]. To make a definite calculation, we need a spin-2 polarization 
vector. An appropriate choice can be found by introducing the massless vectors 1 
and 2 as above and writing 

^mn ^ (1|^-|2] (1|7'^|2] . (29) 

This effectively decays the masive spin-2 particle into two massless spinors. This 
method was introduced in |24j to compute the relevant amplitudes for the emission 
of massive gravitons at high-energy colliders. 



With this prescription, we generate the g-g+ antenna using the operator 

1 
4 



= ^tr[7"^(a-F)7"(a-F)](l|7^^|2](l|7n|2] (30) 
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+ a 



Figure 4: Feynman diagrams for the computation of the qg — )• qgg spUtting functions. 



The matrix element of this operator that creates the 2-parton dipole is 

{g^g^\O\0) = {lAf[2B]\ 



(31) 



To obtain the correct initial polarizations, we take 1 = B, 2 = A as before. The 
matrix elements to the possible 3-parton final states are 



AiO ^ g+g+g+: 
AiO^g+g+g^] 

A{0 g+g-g+: 

A{0 ^ g-g+g+] 







(1&)^[12]' 
{ab) {ac) {cb) 

{ab){ac) {cb) 
(la)^[12]^ 



(32) 



{ab){ac){cb) 

and the conjugates with 1^2 for the other four combinations. Applying the 



reductions (26), (27), we find the results given in the second line of the table. 



The nonzero matrix elements of this operator to qqg final states are 

(Ic)^--^ 



A{0 q+q-g^ 
A{0 q^q+g^ 



ac 



\ac\ 



(33) 



The same reduction process gives the results in the fourth line of the table. 



Spin- 1 and spin-| cases 



The cases of quark-gluon antennae can be treated in the same way. There is one 
additional subtlety. In QCD, quarks are color triplets and gluons are color octets. 
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so a quark-gluon operator carries net color. This means that the matrix element for 
gluon emission from a quark-gluon operator is not gauge-invariant unless we allow the 
gluon also to be emitted from the initial state. This makes it unclear how to define 
a quark-gluon antenna. 

We resolve this problem with the following prescription: We consider the quarks 
to be color octet particles like the gluons. Then, as in the previous sections, we 
extract the color-ordered contribution corresponding to emission from the antenna. 
In the limit of large N^, the various antennae in a process radiate independently. The 
diagrams contributing to a quark-gluon antenna in this prescription are shown in 
Fig. |4j The third diagram, with an intermediate quark line, does not appear in QCD. 
However, it does nicely provide the missing piece that makes this sum of diagrams 
gauge-invariant without radiation from the initial state. 

This solution is the same as that found in the earlier work of Gehrmann-De Ridder, 
Gehrmann, and Glover [TB]. Those authors computed the quark-gluon antennae by 
factorizing the amplitudes for the decay of a neutralino into a massless gluino plus 
gg or qq. In their calculation, the off-shell color octet fermion is the gluino. 

With this understanding, we proceed as in the previous Section. We can generate 
the q-g- antenna using the operator qi{(T ■ F). The polarization spinor can be built 
by introducing massless spinors 1 and 2 as above and taking |2) to be this spinor. 
Then 

= -tqr.ia-F)\2) . (34) 
The matrix element of this operator that creates the 2-parton dipole is 

{q^g^\O\0) = {AB)[B2] . (35) 

To obtain the correct initial polarizations, we take 1 = B, 2 = A. 
The matrix elements to the possible 3-parton final states are 

A{0 ^ q-g+g+) = 



A{0 q^g^g^ 
A{0 ^ q^g+g^ 



{acy{2c) 
{ab){ac){cb) 

{abf{2b) 
(ab) {ac) (cb) 



A(0 ^ t-S-S-) = '-f^^ . (36) 

[ab\ [ac\ [cb\ 



Applying the reductions (26), (27), we find the results given in the seventh line of 
the table. 
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The nonzero matrix elements of this operator to qqq final states are 

. , N (etc) (2c) 

(cb) 

(ab)(2b) , , 

^(0->g_g+g_) = -^^^ . (37) 

The same reduction process gives the results in the ninth line of the table. 

We generate the antenna using the spin-| operator 5^7^(0" ■ F\ This is es- 
sentially the supersymmetry current of the system of gluons and color octet fermions. 
The polarization spinor can be built by introducing massless spinors 1 and 2 as above: 

= zg^7-(a-F)2](l|7^|2] . (38) 

The matrix element of this operator that creates the 2-parton dipole is 

(g_(7+|O|0) = (lA)[2i?]\ (39) 

To obtain the correct initial polarizations, we again take \ = B ^2 = A. 
The matrix elements to the possible 3-parton final states are 



A{0 q-g+g+] 
A{0 ^ q-g-g+] 
A(p q^g+g^] 



{ab) {ac) {cb) 
{ab)[2bf{12) 

[ab] [ac] [cb] 
{ac)[2c]'{12) 



[ab] [ac] [cb] 

A{O^q^g^g^)=0 . (40) 



Applying the reductions (26), (27), we find the results given in the eighth line of the 
table. 

The nonzero matrix elements of this operator to qqq final states are 

{la)[2b]' 



A{0 q-q-q^ 



[cb] 



AO^q^-q^q^)^-^^. (41) 
The same reduction process gives the results in the tenth line of the table. 
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6 Initial-state showers 



The Feynman diagram computations that we have done to find the antenna spht- 
ting functions for FF sphttings can also be apphed, by crossing, to IF and II split- 
tings. The expressions in Table [T] are given in terms of invariant quantities that are 
unchanged under crossing. Thus, we can use the expressions in this table directly in 
other channels. At worst, a change of the overall sign is required in some cases. In 
this section, we will clarify this statement by analyzing the kinematics of IF and II 
splittings in the same variables as those used in Section 2 for FF splittings. In all 
cases, the kinematics is done for all massless partons only. The kinematic discussion 
in this section is similar to that presented in [23j. 

To begin, we will formalize some of the results quoted in Section 2 for the FF 
region. The cross section for a process X — )■ ach is 

a{X -> ach) = :^^T-^ / dZadz^\M{X -> ach)\^ , (42) 

where $x is the flux factor. Polarization and color indices have been suppressed. The 
left-hand side has been integrated over the orientation of the final state system but 
is otherwise exact. To write an expression involving the antenna splitting function, 
we approximate 

M{X^ ach) ^M{X^AB).gT. ^ , (43) 

where O is the operator used in Sections 3-5 to represent the state AB. The factor 
gT is the QCD coupling and color matrix; after squaring and summing over colors. 



this becomes ATragNc- The splitting function is defined by (13), 

M{0^ach) " 



S{Za,Zc,Zb) = SAB 

Then 



M{0 AB) 



(44) 



a{X ach) ^ a{X -> AB) ■ — — - [ dzadzbS{za, z^, Zb) ■ (45) 

47r J 



It is important to note that, in this formula or in (43), the vectors kA and /c^ are 
introduced as part of the approximation. They can be defined in any way that is 
consistent with the requirements that kA and ks are lightlike, kA + ks = Q, and kA 
and ks become parallel to ka and kb, respectively, in the soft and collinear limits. 

The logic of this derivation extends straightforwardly to the IF and II regions. 
The major change is that, in these cases, we need to introduce initial hadrons from 
which the initial partons are extracted. 
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Consider first the IF case. The cross section for a proton of momentum P to 
scatter from a color-singlet system X transferring momentum Q to create a 2-parton 
system ch is 

aipX -^cb)= I dxafixa) ^— I dcosd^ \MiaX ^ c6)|% (46) 

J ^ax 167r J 

where cos ^* is the scattering angle in the ch center of mass system. We will approxi- 



mate this formula using the expression analogous to (43) 



M{aX ^ ch) ^M{AX^B).gT. ^[""f^^"^] ■ (47) 



Then the splitting function is defined by the same expression S as in (44), but now 
analytically continued into the new kinematic region. If a fermion line is crossed from 
the final to the initial state, an extra factor (-1) should be included. In addition, sab 



in (44) is negative in this region, giving an extra minus sign. 



The decomposition of the amplitude is illustrated in Figjsta). The kinematics 
can be described by variables yij and Zi obeying the relations ([l| to However, the 
vectors fc^, ka now have negative timelike component, and the vector Q = kA + ks = 
ka + kf, + kc is spacelike, = sab < 0. The phase space for this region covers the 
quadrilateral shown in Fig. |5](b). The region of integration is infinite, since Za can 
become very large, but the integral is cut off at large Za by the parton distribution 
function. The line Za > 1, Zf, = 1 corresponds to the region of initial state radiation, 
c parallel to a. The line Za = I, < Zb < I corresponds to the region of final 
state radiation, c parallel to h. The line Za + Zb = 1 corresponds to h parallel to a, 
that is, h as initial state radiation from the primary a. An antenna shower should 
give an accurate description of the dynamics in the two regions \yac\ < \ybc\ < 1, 
lUbcl < IVacl < IVabl that are shaded in the figure. The new constraint \ybc\ < 1 is just 
\sbc\ < \Q^\, which is stronger than the constraint that this invariant is less than \sab\- 



To decompose (46 ) into an appropriate form, we choose pa and and then change 
variables. Let pA be chosen in the direction of pa, so that pa = ZaPA, Za > I. Then 
Pb = Q — Pa- We have 

Pa = XaP , Pa = XaP , so Xa = ZaXA , (48) 

with xa having the definite value xa = —Q^/2P ■ Q associated with scattering a 
massless particle from a local current. For the reaction aQ he, s + t + u = Q'^, so 
t + u = Q'^-s = Q'^Za. Then 

t = Q\l - Zb) = ^Q^Za{l - cos 9,) (49) 
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Figure 5: (a) Kinematics of 2 — )• 3 splitting in the initial-final (IF) case, (b) Phase space 
for 2 — )• 3 splitting in the IF case. The eight regions corresponding to different orderings 
of \sab\, \sac\, Sbc, \Q'^\ are shown. The region that should be well described by an antenna 
splitting AB — >• acb is shaded. 



Using these formulae, we can change variables from (xa,cos6'*) to {za,Zb). The Jaco- 
bian of this transformation is 

_ d{xa,cos 9^) _ 2xa .^„s 

d{Za,Zb) Za 

Thus, 

a{pX cb) = J ~^ J '^^b J dXAXAfiZaXA) S{xa + Q'^/2P ■ Q) 

■^hM{aX ^ ch)\\ (51) 



This is an exact rewriting of (46). Now apply the approximation (47) and group 
terms to form 

a{AX ^B) = -^27r5(Q2 + xa2P ■ Q)\M{AX ^ B)\^ . (52) 

Then 

a{pX Cb) ^ I ^ I I (^^AfiZaXA)(T{AX B) ■ ^^^^S{Za,Zc,Zb) . (53) 

As an example, consider using this formula to describe initial-state gluon radiation 
in deep inelastic scattering from a quark. The total gluon emission is given by the 
sum of the two spin-dependent splitting functions in the fifth line of Table [T} equal to 

^5 = _^i±ff, (54) 

yacVch 
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The extra minus sign comes from the sign of sab in (44). In the region of initial state 



radiation, = l/w, ^ 1, Vac = -(1 - ^fe), Vbc = (1 - V^)- Then, setting 



1 

dzb- = log ^ , (55) 

1- Zb /i^ 



we obtain 



a{pX ->cb)^ / dxA / — / — criAX B) ■ — — -log — , 56 

J J W W 4:7T [1 — W] fi^ 



which is correct. 



This is an appropriate point to discuss again the signs of the expressions in Ta- 
ble [T} The antenna splitting functions are probabilities; thus, they should be positive. 
However, we define the splitting functions in the IF and II regions as analytic contin- 
uations of the values in the FF region, so their positivity must be checked explicitly. 

As we move from the FF region to the IF region with A and a in the initial state, ybc 
becomes negative while all other yab, Vac, Za, Zb remain positive. The factor Zc can be 



negative, but Zc does not appear in the Table. With the minus sign from sab in (44), 
the denominator of S{za, Zb, Zc) is positive, and so we need only check the numerator 
functions in given in the Table. The numerator functions for gg — )■ ggg, qq — )■ qgq, 
and qg — > qqq remain positive, while the numerator functions for gg — )■ qqg become 
negative. In this last fermion not present in the 2-parton system is crossed 

from the final to the initial state, so we must supply an extra factor (—1). Then all 
of the expressions are positive, as required. However, if we then cross from the region 
26 > to the region Zb < 0, one qg — )■ qgg and one qg — )■ qqq amphtude changes 
sign. This sign change is unphysical; presumably, it is due to the simple method of 



reconstruction in (26) and (28). We recommend setting these two amplitudes to zero 
for Zb < 0. The region 2;;, < is outside the shaded region in Fig. |5]^b) where we will 
generally use the parton shower approximation, so most likely this difficulty is not 
important in practice. 

Similarly, for the FI region where b and B and taken to be in the initial state, 
i/ac < 0. Then the numerators that go negative as we cross into the region are those 
in the qg — )■ qqq cases where a fermion is crossed into the initial state. Now there are 
four amplitudes, one each in the qg — t- qgg cases and both of those in q-g+ — qqq, 
that become negative when Za < 0. Again, we recommend that these amplitudes be 
set to zero in this region of unphysical behavior. 

In the II region, both i/ac and i/bc are negative. The denominator of S{za, Zb, Zc) 
is positive. The numerator terms that are negative because of the sign changes are 
compensated by minus signs from crossing. There are no unphysical sign changes. 
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Figure 6: (a) Kinematics of 2 — )• 3 splitting in the initial state (II) case, (b) Phase space 
for 2 — )• 3 splitting in the II case. The six regions corresponding to different orderings of 
I Sac I; \sbc\, IQ^I are shown. The region that should be well described by an antenna splitting 
AB — >■ acb is shaded. 



The correct result is always obtained by taking the absolute value of the numerator 
expression from Table [Tj 

We now discuss the kinematics of the II case. We begin from the formula for two 
protons of momentum P4, Pg to produce a color-singlet system of momentum Q plus 
a massless parton c, 



a{pp cX) = [ dxa [ dxbf{xa)f{xb) 7^—7^ [ d COS 9^—^ \M{ab ^ cX)\ 

J J 2s ah IbTC J . Sab 



11/-,, 

(57) 

where cos 9^ and are the scattering angle and the momentum in the cX center of 
mass frame. 



The decomposition of the amplitude is illustrated in Fig. |6[a). The kinematics 
can again be described by variables yij and Zi obeying the relations ([T]) to (|4]). Now 
the vectors kA-, ka ks, kf^ have negative timelike component, and the vector Q = 
kA + ks = ka + kb + kc is also negative timelike, with > 0. The phase space 
for this region covers the quadrant shown in Fig. |5](b), with Za, Zb > 1. Again, the 
region of integration is infinite, but the integral is cut off by the behavior of the 
parton distribution functions. The hne Za > I, Zb = 1 corresponds to the region of 
initial state radiation with c parallel to a. The line Za = 1, zi, > 1 corresponds to 
the region of initial state radiation with c parallel to b. An antenna shower should 
give an accurate description of the dynamics in the two regions \yac\ < \ybc\ < 1, 
\ybc\ < \yac\ < 1 that are shaded in the figure. Again, the limit 1 here corresponds to 
constraints |sac|; l^hd < \Q'^\, which are stronger than the constraints that these two 
invariants are less than \sab\- 
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In the ab — )■ cX process, the system X must recoil with some nonzero transverse 
momentum. Thus, it is not possible to choose fc^ and fc^ to be parallel to ka, kh. 
The invariants for the ab — )■ cX scattering process satisfy s + t + u = Q^. Since 
t = — Zb), u = Q'^{1 — Za), this means that s = Q^{za + Zh — 1). Alternatively, 
s = XaXb ■ IP A ■ Pb- We would like to choose the longitudinal fractions of A and B, 
Xa and xs, to satisfy the relation 

XAXB ■ 2Pa ■Pb = Q^ . (58) 
To make this possible, we must write 

Xa = ZaXAC , Xb = ZbXsC , (59) 

with [25] 

= ^ ~ ^ (60) 

ZaZb 

The function C{za, Zb) approaches 1 when either Za or Zb goes to 1; that is C ~ 1 in 
both collinear regions. 

Also, t + u = Q"^ {2 — Za — Zb) = Q'^Zc, so 

t = Q\l- Zb) = \q^z,{\ - COS d,) . (61) 



We can now use (59) and (61) to change variables from (x^, Xf,, cos6'*) to (x^, 2;^, z;,), 
holding Xb fixed at the value xb = Q"^ / xa2Pa ■ Pb- The Jacobian of this transforma- 
tion is 

d{xa,Xb,cose^) 2xb Q"^ Js^b 
J=— ^ ^ = = xb — . (62) 

d[XA,Za,Zb) Zc Sab P* 



Then 



a{pp cX) = J J dxAdxB f{zaXAC)f{zbXBC) xbS{xb - Q'^/xa2Pa ■ Pb) 



■ — — \M{aX^cb)\\ (63) 
Sab ott 



This is an exact rewriting of (57). Now apply the approximation analogous to (43) 
or (47) and group terms to form 

a{AB ^X) = -^2n6{Q^ - xaXb2Pa ■ Pb)\M{AX B)\\ (64) 
2sab 

This gives, finally, 

a{pp cX) ^ J J dXAdXB f{ZaXAC)f{ZbXBC)a{AB X)-'^^S{Za, Zc, Zb) 

(65) 
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To test this formula, consider the case of qq annihilation with the emission of a 
gluon collinear with the quark a. The sum of spin-dependent splitting functions for 
this case is again (54). In the collinear region of interest, = z;, ~ 1. Repeating 
the step that led to ( 56 ) , we find 

dw 



a{pp -> cX) ^ / dxAdxB [ —f{—)f{xB) (y{AB X) ■ ^ x log ^ 

J J w w An [1 — w) 

which is the correct limit. 



(66) 



7 Comparison to previous results 

In the Introduction, we made reference to a number of previous definitions of the 
antenna splitting functions. We noted that these definitions agree, as they must, in 
the singular soft and collinear limits. However, these prescriptions differ widely away 
from the boundaries of phase space. In this section, we will compare our prescription 
to those of ARIADNE id and Gehrmann-De Ridder, et al. [IS]. 



We will make this comparison over the natural phase space discussed in the pre- 
vious section-the entire {za, z^) plane above the line Za + Zb = 1. In order to describe 
antenna showers for initial- as well as final-state emissions, the splitting functions 
should extend into the region Za,Zb > 1. Depending on the details of how the shower 
is constructed, their use might be restricted to a polygon around Za = Zb = I, or the 
expressions might be used for arbitrarily large values of Za and Zb- 

We note again that the IF regions include the lines Za = and Zb = 0. Expressions 
for the splitting functions that are well-behaved near Za = Zb = 1 can possibly have 
a singularity on this line, though such a singularity in the middle of the phase space 
would be unphysical. We used this criterion in Section 4 to exclude factors of 1/za 



and 1/zb from appearing in (28). The antenna functions of Duhr and Maltoni 
are typically singular along this line and so cannot be used in parton shower models 
in all regions. 

The ARIADNE and Gehrmann-De Ridder antenna functions give expressions 
summed over final polarizations. To compare our splitting functions to these, we 
must sum over a row in Table [1} Our summed expressions are independent of the ini- 
tial polarization in the soft and collinear limits, but they depend on the polarizations 
of A and B in the interior of the [za, Zb) space. The comparison to our expressions 
thus also reveals where this dependence on polarization is an important effect. 

The first antenna splitting functions were put forward by the ARIADNE group [6] . 
Their approach started from the spin-averaged cross section for the simple splitting 
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process qq — )■ qgq in e^e~ annihilation. They then guessed the expressions for the 
~^ (199 S'l^d gg ggg sphttings, so that these would have a similar form to the 
qq — i- qgq case, 

S = ° ^ ' , (67) 

yacVbc 

where Ua, Ub = 2 for emission from a quark and 3 for emission from a gluon. 

Our philosophy, explained in Section 2, is that each individual antenna should 
reproduce the coUinear limit predicted by QCD. These expressions are symmetric 



under interchange of identical particles, while (67) does not have this property, so we 



would obtain the complete splitting function by symmetrizing (67). This gives 
qq antenna: S 



4 + 4 






VacVhc 






4 + zl 


+ 


zl + zl 


VacVbc 


VabVbc 


4 + 4 


+ 


zl + 4 


VacVbc 


VabVbc 



gg antenna: S = — 1 — + 



zl + zl 



VabVac 

qg antenna: S = ~° ' + ~" ' . (68) 

VacVbc VabVbc 

The summed terms are each positive in the FF kinematic region. To obtain the 
ARIADNE splitting functions in the other regions, we analytically continue these 
formulae into the regions where Za or Zb is greater than 1. 

The analytic continuation of the ARIADNE and, below, the Gehrmann-de Ritter 
results brings in the issue of the positivity of these expressions, similar to the pos- 
itivity issue for our splitting functions discussed in Section 6. For the ARIADNE 
and Gehrmann-De Ridder antenna functions, the expressions given are summed over 
spins, and the individual pieces are not independent of one another. So, if they be- 
come negative, that is a problem for the complete, spin-summed, expression. For 
the Gehrmann-de Ridder functions, it can be seen that this happens only the regions 
Za < and Zb < 0, so this is not a serious problem. However, the ARIADNE function 
involve zl, which is negative in the whole region Za + Zb > 2. This problem cannot 
be resolved by replacing Zc with \zc\, since this leads to expressions that do not agree 
with the Altarelli-Parisi factorization along the lines separating the IF regions from 
the II region. Fortunately, the ARIADNE functions do not become actually become 
negative until Za or Zb becomes very large {za or Zb ~ 12). However, the idea that 
the ARIADNE functions are sums of positive and negative terms in the initial-state 
regions goes against the intuition used to propose these expressions. 

We are now in a position to compare the ARIADNE function to our proposal. 
For the qq antenna, the expression above coincides with the sum of row 5 of Table 
[TJ For the gg and gq cases, the ratio of the above ARIADNE functions to those 
defined in Table [T] are illustrated in Figs. [7| [8| and [9j The notation in the figures is 
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2 

210123 210123 

Figure 7: Visualization of the ratio of the ARIADNE antenna function to our antenna 
functions for the processes gg — t- ggg. The figures on the left and right are the comparison 
of the ARIADNE antenna function to our spin-summed antenna functions from row 1 and 
row 2 in Table [T| respectively. The boundaries of phase space for the different kinematic 
regions are marked in blue. The contours are plotted at ratios of 1.2, 1.5, 2.0, 3.0, and 5.0, 
with + indicating a region in which the ratio is greater than 1. 



the following: Each figure represents the ratio of the ARIADNE splitting function 
to our results for a specific initial set of polarized partons, summed over final state 
polarizations. The ratio goes to 1 on the lines Za = 1 and Zb = 1, which correspond 
to the collinear limits. Away from these lines, the contours on which the ratios are 
1.2, 1.5, 2.0, 3.0, and 5.0 (toward the + symbol), and the inverses of these numbers 
(toward the — symbol) are shown. The qg antenna function are asymmetric between 
partons a and b. The IF region in the lower right is that in which the quark is in the 
intial state and the gluon is in the final state. The IF region in the upper left is that 
in which the gluon is in the initial state and the quark remains in the final state. 



The ARIADNE authors gave a different interpretation to the formulae (68). They 
took the philosophy that the collinear limit need not result from a single antenna but 
rather should be the result of summing over the possible antennae that would lead 
to a specific final state. A three gluon final state could result from any pair of the 
gluons radiating the third and so should be the sum of three antennae. Then the 



second line of (68) would be interpretated as the sum over these three antennae. This 
is a reasonable point of view for the FF kinematics considered in [6] . However, in the 
IF and II regions, at least one of the Zi will be negative and so some of the terms 
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Za 



Figure 8: Visualization of the ratio of the ARIADNE antenna function to our antenna 
function for the process q-q_ — )• qgq. Our antenna function for the process q-qj^ — )• qgq 
coincides with the ARIADNE result and so is not included. The notation is as in Fig. [7| 




210123 210123 

Za Za 

Figure 9: Visualization of the ratio of the ARIADNE antenna function to our antenna 
functions for the processes qg — )• qgg. The figures on the left and right are the comparison 
of the ARIADNE antenna function to our spin-summed antenna functions from row 7 and 
row 8 in Table [T| respectively. The notation is as in Fig. [7j 
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in the last two lines of (68) will become negative. Such terms cannot be interpreted 
as independent radiators, each emittting a gluon with positive probability. It is 



tempting to revise the formula in (68) by taking the absolute values of the negative 
terms. However, one can readily check that no such prescription gives the correct 
Altarelli-Parisi limit along the lines Za = 1 and zi, = 1 ai the boundaries of the IF and 
II regions. Thus, we believe, the ARIADNE formulae can be used in the IF and II 



regions only by using the formulae (68) as written and accepting that some negative 
signs will appear 



Gehrmann-De Ridder, Gehrmann and Glover [TB] studied 2 — )■ 3 splitting from 
Feynman diagrams to develop an antenna subtraction program for NNLO calcula- 
tions. In doing so, they were able to extract unpolarized antenna functions for the 
processes gg — )■ ggg, qg — >■ qgg and qg — )■ qqq. To calculate the gluon-gluon antenna 
function, they used the effective Higgs coupling to gluons 



C = -^hF'"'F. 



4 



(69) 



This is essentially the same procedure that we used in Section 3, and it yields the 
same result as the sum of row 1 in Table [Tj In our language, their antenna function 
for the gluon-gluon dipole is [27] 



S 



vie + vie + Vlh + vWhc + vlhVlc + vlhVlc 



+ 4 . 



(70) 



VahVacVhc 

The comparison of this antenna function to the sum of row 2 of Table [T] is illustrated 
in Fig. [To 



This splitting function for gg — )■ ggg is, however, not precisely the form of the 
splitting function that is used in the VINCIA parton shower [8]. They use the 'global' 
form of the Gehrmann-De Ridder antenna function, which in our language is 



S 



'^ylb + yabvlc + vlbylc , 8 

VabVacVbc 3 



(71) 



To implement this antenna function, a similar procedure is used as with the ARI- 
ADNE antenna functions. That is, emissions from overlapping antenna are summed. 
When the three antennae contributing to gg — )■ ggg are summed together, one recov- 



ers the result (70). This prescription works well in the FF kinematics. However, as 



in the ARIADNE case, it might require negative contributions in splitting functions 
for some antennae in the IF and II kinematics. 

To construct the antenna functions involving quarks, Gerhmann-De Ridder, et ai, 
calculated the decay of a neutralino x to a gluon and a gluino ip through the effective 
operator 



C = ti^^a^^xF^, + h.c 



(72) 
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Figure 10: Visualization of the ratio of the Gehrmann-De Ridder antenna function to 
our antenna function for the process g-g-\- — )• ggg. The antenna function for the process 
g+g+ —7- ggg coincides with the Gehrmann-De Ridder result and so is not included. The 
notation is as in Fig. [7j 



In principle, our results should agree for case of a spin \ initial state. However, our 



choices (26) and (28) for handling ambiguous momentum products, produce some 



differences. In our language, their antenna functions involving quarks are 



yabVacybc 



qg qqq: S 



iVac + VabfyacVab - '^vlcVlb 



yabVacVbc 



+ Vab + Va 



(73) 



The comparison to our antenna functions is illustrated in Figs. 11 and 12 For 



qg — )■ qgg, our result for the spin \ case is indeed very close to the above expression 



in the FF region. For qg — )• ggg, our prescription (28) gives us an extra factor of Za 
0. 



near z, 



a 



In summary, we have shown that the antenna splitting functions represented by 
^ and Tabled give a physically sensible prescription for the construction of antenna 
showers. These splitting functions can be used with the formulae (45), (53), (65) to 



generate antenna splittings in all three relevant kinematic regions. We hope that this 
formalism will provide a firm foundation for the construction of new parton showers 
based on the antenna concept. 
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Figure 11: Visualization of the ratio of the Gehrmann-De Ridder antenna functions to 
our antenna functions for the processes qg — )• qgg. The figures on the left and right are 
the comparison of the Gehrmann-De Ridder antenna function to our spin-summed antenna 
functions from row 7 and row 8 in Table [T| respectively. The notation is as in Fig. [7] 




Figure 12: Visualization of the ratio of the Gehrmann-De Ridder antenna functions to 
our antenna functions for the processes qg — >• qqq. The figures on the left and right are 
the comparison of the Gehrmann-De Ridder antenna function to our spin-summed antenna 
functions from row 9 and row 10 in Table [T| respectively. The notation is as in Fig. [7| 



27 



ACKNOWLEDGEMENTS 



We thank Darren Forde, Tanju Gleisberg, Peter Skands, and Jan Winter for in- 
structive discussions and Claude Duhr and Fabio Maltoni for a useful correspondence. 
This work was aided by our participation in the Northwest Terascale Workshop on 
Parton Showers and Event Structure at the LHC at the University of Oregon. We 
thank the participants and, especially, the organizer, Davison Soper. The work was 
supported by the US Department of Energy under contract DE-AC02-76SF00515. 



References 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 

[7 



[9 
[10 

ill 



T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) |arXiv:hep- 
PV0603175j . 



G. Corcella et ai, JHEP 0101, 010 (2001) [arXiv:hep-ph/0011363| , [arXiv:hep- 



ph/0210213 



G. Altarelh and G. Parisi, Nucl. Phys. B 126, 298 (1977). 

Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977) [Zh. Eksp. Teor. Fiz. 73, 1216 
(1977)]. 

G. Marchesini and B. R. Webber, Nucl. Phys. B 238, 1 (1984). 

G. Gustafson and U. Pettersson, Nucl. Phys. B 306, 746 (1988); B. Andersson, 
G. Gustafson and L. Lonnblad, Nucl. Phys. B 339, 393 (1990). 

U. Pettersson, Lund preprint LU-TP-88-5 (1988); L. Lonnblad, Comput. Phys. 
Commun. 71, 15 (1992). 

W. T. Giele, D. A. Kosower and P. Z. Skands, Phys. Rev. D 78, 014026 (2008) 
|arXiv:0707.3652 [hep-ph]]. 



J. C. Winter and F. Krauss, JHEP 0807, 040 (2008) |arXiv:0712.39l3| [hep-ph]]. 
M. E. Peskin, in preparation. 



S. Catani and M. H. Seymour, Phys. Lett. B 378, 287 (1996), |arXiv:hep- 



ph/96 02277]. Nucl. P hys. B 485, 291 (1997) [Erratum-ibid. B 510, 503 (1998)] 
" |a5Civ:hep-ph/ 9605323| . 



[12] D. A. Kosower, Phys. Rev. D 57, 5410 (1998) |arXi v: hep-ph/ 971 02 13| . 



28 



[13] J. M. Campb ell and E. W. N. Glover, Nucl. Phys. B 527, 264 (1998) |arXiv:hep- 
lph/9710255]. 

[14] S. Catani and M. Grazzini, Phys. Lett. B 446, 143 (1999) [arXiv:hep- 
pli/981038 9]. 

[15] C. Duhr and F. Maltoni, JHEP 0811, 002 (2008) |a rXiv:0808.3'3l9| [hep-pli]]. 

[16] A. Gehrmann-De Ridder, T. Gehrmann and E. W. N. Glover, Phys. Lett. B 
612, 49 (2005) |arXiv:hep-ph/0502110], Phys. Lett. B 612, 36 (2005) |arXiv:hep- 
' ph/050129 1], JHEP 0509, 056 (2005) [arXiv:hep-ph/0505111] . 

[17] In writing a parton shower for e~^e~ — )■ qq, one might pick up an obvious 1/Nc 
correction by writing the prefactor as | instead of 3 for the first gluon emission 
from the quark line. In higher orders, the analogous treatment of radiation from 
quark lines is not so innocent and can potentially lead to incorrect accounting of 
color coherence. 

[18] M. L. Mangano and S. J. Parke, Phys. Kept. 200, 301 (1991) |arXiv:hep- 
th/050922J3j. 

[19] L. J. Dixon, in QCD and Beyond: TASI 1995, D. E. Soper, ed. (World Scientific, 
1996) arXiv:hep-ph/9601359 . 

[20] L. J. Dixon, E. W. N. Glover and V. V. Khoze, JHEP 0412, 015 (2004) 
[arXiv:hep-th/0411092] . 

[21] R. Kleiss and W. J. Stirling, Nucl. Phys. B 262, 235 (1985). 

[22] D. A. Kosower, Phys. Rev. D 67, 116003 (2003) )arXiv:hep-ph/0212097|. 

[23] A. Daleo, T. Gehrmann and D. Maitre, JHEP 0704, 016 (2007) |arXiv:hep- 
ph/0612257|. 

[24] M. P. Le and M. E. Peskin, in preparation. 

[25] Compare eq. (5.15) of [23j. 

[26] The problematical terms appear to have been simply omitted in the parton 
shower model of [U]. We thank Jan Winter for an extensive discussion of this 
point. 

[27] This expression actually differs from the one given in [16J by a factor of 1/3, 
which comes from allowing any gluon to become coUinear with any other gluon. 
Since we have identified the radiated gluon, we remove this factor. 



29 



